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Abstract 



In this paper, we analyse the optical absorption in porous 
silicon. This is the first attempt to explicitly demonstrate that it is not 
possible to extract the band gap of such low dimensional nanostructures like 
porous silicon from a Tauc plot of V otfiuj vs fiuj. So we model the absorption 
process assuming that porous silicon is a pseudo ID material system having 
a distribution of band gaps. We show that in order to explain the absorption 
we specifically need to invoke - (a) k is not conserved in optical transitions, 
(b) the oscillator strength of these transitions depends on the size of the 
nanostructure in which absorption takes place and (c) the distribution of band 
gaps significantly influences the optical absorption. A natural explanation of 
the temperature dependence of absorption in porous silicon also follows from 
our model. 
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I. INTRODUCTION 



The discovery of efficient photoluminescence ^^'^^ from porous silicon (PS) has 
attracted the attention of many researchers. The photoluminesecence peak occurs at an 
energy greater than the band gap of c-sihcon and can be tuned through the visible 
spectrum by changing the preparation conditions. Canham proposed that PS is a 
nanostructured material and the band gap of these nanostructures is enhanced due to 
quantum size effects. He argued that this would account for the fact that the luminescence 
energy is greater than the band gap of c-silicon. Evidence for the increase in band gap was 
first reported by Lehmann - from optical absorption measurements Transmission 
electron microscopy measurements reveal that the diameter of the columnar 
nanostructures that make up PS depends on preparation conditions and can be as small as 
2-4 nm PS can hence be thought of as an assembly of pseudo ID quantum wires. 
The progress in this area has been the subject of many recent reviews In this paper 

we confine ourselves to the understanding of the optical absorption process in PS. 

There have been many attempts to obtain the band gap of PS from optical 
absorption measurements In these measurements, the authors assume that the 

absorption coefficient {a) of PS satisfies the same relation as in 3D bulk c-silicon viz. 

= A{huj - Eg) (1) 

There the Tauc plot of \J olTho vs hjj is a straight line and the intercept on the 
energy axis gives the band gap Eg. The use of this equation for a 1 dimensional system like 
PS suffers from the following criticisms . 

1) In general ,the absorption co-efficient depends on the joint density of states 
(JDOS) of the material. Eq.l is valid only when the density of states g(E) varies ^^'^^ with 
energy near the band edges as g(E) = M-\/E', where E is measured from the band edge and 
M is a constant. If we treat PS as a ID material ( upto a first approximation assuming 
parabolic band structure ) it would be more appropriate to write the density of state g(E) 
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= ( N is a constatnt ) than the common practice of using 3D density of state. Clearly 
the JDOS in PS is expected to be different from that of bulk c-silicon. For a ID system 
having a single direct band gap Eg, it is easy to show that the JDOS is proportional to 

, In such plot of a{E) vs E should peak at Eg. On the otherhand, for an 

indirect band gap material ( assuming that the momentum k is not conserved in this 
pseudo ID system in an optical transition ), it is easy to show that the JDOS is a constant 
and independent of incident energy ( following reference (14) ) . Therefore , in either case , 
\/ahio is not linearly related with hw at all. Clearly the use of eq. 1 to extract the band 
gap is not justified for PS . 

2) The nano structures that make up PS have a distrihutiorS^'^'^^ of diameters 
(d). If the band gap is related to quantum size effects, then it is clear that the energy 
upshift ^E ( which can be written as AE = C/d^ ) is different for each of the 
nanostructures that make up PS.( The value of X in the simplest approximation is 2 , 
further details of this are discussed in section III ). In PS, we are actually deahng with a 
heterogeneous system with a distribution of band gaps. Hence it may not be meaningful to 
visualize PS as having a single band gap and that this can be extracted from equation 1. 

At this junction, we would also like to point out that the above objections 
against the common practice of using ^^^^ the eqn. 1 to get the band gap are equally vaild 
in case of other low dimensional systems like quantum dots etc. 

In an attempt to address these questions in case of PS, we have investigated the 
optical absorption in PS specifically using the ID density of states . In this paper we 
clarify the absorption process in PS assuming that PS has a distribution of band gaps. 
Section II contains the experimental results which will be used for comparing our 
simulations. Here we report a non — destructive way of measuring the porosity of PS thin 
films using transmission measurements and effective mass theory. Section III outlines the 
model used for calculating optical absorption and section IV the results of our simulations 
. On the basis of the simulations we will further consolidate our viewpoint that the Tauc 
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plot can not be used to determine the band gap of such low dimensional systems . We also 
show that we can easily explain the temperature dependence of the absorption in PS using 
our model .Finally we summarise our conclusions in section V . 



II. EXPERIMENTAL RESULTS 

We have done some experiments to facilitates the comparison of our simulations 
with real experimental data . 

PS is made by electrochemical etching of p type c-si in a 1:1 HF : Ethanol 
solution at a current density between 10 - 50 mA/cm^. Free standing samples are lifted off 
the substrate by increasing the current density to > 200 mA/cm^. The samples are then 
washed in water to get rid of HF , rinsed in n-pentane and finally dried in air on a glass 
substrate for transmission measurements . 

Transmission measurements are carried out using a CVI 240 Digikrom 
monochromator with a tungsten lamp as a source, silicon detector and SR 530 lock -in 
amphfier .For some samples we also use a Gary 1756 spectrometer. We have checked that 
internal multiple scattering does not dominate the absorption in our samples. This is done 
by making two PS samples of two different thickness (5// and 10// respectively) and 
ensuring that the transmission curves are the same for both samples at low absorption ^^^^ . 

The absorption coefficient a (A) as a function of wavelength A is obtained from 
the normalised transmittance T(A) using ^^^^ 

^ ' 1 - R^exp{-2a{X){l - P)t) ^' 
The reflectivity (R) is obtained in the low a region using the relation R = 

(1-T)/(1-|-T). P is the porosity of the layer and t is its thickness. 

The porosity of PS is usually determined by gravimetric measurements which 

also destroy the sample. We now show that optical transmission experiment allows us to 

measure porosity in a nondestructive fashion . 
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The normalized transmittance ( in eq.2) for o; = reduces to 



If n is the real part of the refractive index , then we can also write for the 
reflectivity ( for o; = ) as 

R = 'I (4) 
It follows from eq.4 that the refractive index (ups) of PS is 

We model PS as consisting of two media - air ( having a refractive index ( Uair 
= 1 )) and crystalline silicon ( nsi = 3.44 ) . It was shown earlier ^^^^ that in the 
framework of effective medium approximation one can write , 

{eps - l)/(3ep5) = (1 - P)[{es^ - l)/{2eps + esi)] (6) 
where P is the porosity and the dielectric constant e — Therefore, 



P = 1 - [{2eps + esi){eps - l)/{esi - l)(3ep5)] (7) 

The porosity can now be obtained using equations 4 to 7. 

To check the vahdity of eq.7, we calculate the porosity ( using eq.7 ) and 
compare it with the corresponding values of P obtained gravimetrically. Figure 1 is a plot 
of the porosity determined from the transmission curves and the corresponding 
experimentally determined porosity by gravimetric measurements as reported in the 
hterature*^^'-*^^^ . We see from figure 1 that the points nearly fit a straight hue with slope 
one . We hence conclude that optical transmission measurements can reliably be used to 
obtain the porosity of porous silicon in a non — destructive fashion . 
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Figure 2 is a plot of log(Q;) vs E at 300 K and at 100 K. The absorption at low- 
temperature is reduced with respect to that at room temperature. We see that the two 
curves exhibit a rigid shift in agreement with other results ^^^^ . An explanation of this 
phenomena on the basis of our simulated results is given in section IV . 

We now attempt to understand the optical absorption process in PS with the 
help of some model calculations and simulations, assuming that these ID materials have a 
distribution of band gaps . 

III. A MODEL FOR ABSORPTION IN PS HAVING A DISTRIBUTION OF 

BAND GAPS 

We model the nanostructures that make up PS to be ID parallelopipeds of 
square cross sections of side d and constant length L. Since d is typically 20 - 40 A, the 
band gap of the material is enhanced due to quantum size effects . We assume that the 
effective band gap ( lowest energy gap ) for a particular parallelopiped is the minimum 
separation between the conduction band and the valence band states for quantum numbers 
= 1 and niy = 1. We also assume that the absorption takes place in the bulk. 

The absorption coefficient a{E) for a particular value of incident energy E can 
be written as 

a(E) = / p'(E - EG)F{EG)P{EG)dEG (8) 

JEg<E 

where the integration is done over all the band gaps (Eg) of the system below the incident 
energy E. In equation 8, pj is the joint density of states, F(E(3) the oscillator strength for 
the optical transition and P(E(;) the distribution of band gaps E^ in the material. We now 
discuss eq.8 in some detail. 

We have already pointed out in section I, if g(E) —B/y/E then the JDOS is 
p"^ = A(E — Eg)~^ where A is a constant. The value of W depends on whether PS is a 
direct or an indirect band gap material. For direct band gap material W = 0.5 and for 
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indirect band gap material W = . 

F(Eg) is the oscillator strength governing the absorption. For indirect band gap 
material, the oscillator strength ^^^^ can be enhanced for small d . This happens because 
the overlap between the electron and hole wavefunctions in k space increases as a result of 
quantum confinement and contributes to an increase in the oscillator strength at small d 
(large Eq) and can be written (^s-^^) as F = i/d? where 7 = 5 to 6. This manifests itself as 
the dominant no phonon line which has been seen recently in the luminescence spectra of 
PS 

The third term in eq.8 viz. P{Eg), is the distribution of band gaps in PS. The 
band gap distribution is obtained by assuming a distribution of sizes for d and a relation 
governing the upshift in energy A E with the size d ( due to quantum confinement ). 

In our model we have considered two possible distribution of sizes. One is 
Gaussian and the other is Lognormal^^^^ distribution given by P'^(d) and -P^(d) 
respectively. 



where do is the mean size and the a the standard deviation. 



P'id) = V ea;p[- ^^^q-r°^' ] (10) 



where mo = In(do), and ctl = ln((7), do being the mean size and a the standard deviation. 

Electron microscopy measurements^^'^'^^^ suggest that for p PS, do is around 30 
A and cr ~ 4 A. Figure 3 shows the distribution of sizes for Lognormal and Gaussian 
distributions for do = 30 A and a — A A. 

The energy up shift for the confinement AE can be written as 



AE^Eo-Eg^ C/d^ (11) 
where E^ is the crystalline Si fundamental indirect band gap ~ 1.17 eV and E^ is the 
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increased band gap due to quantum confinement. 

The value of X is 2 in the usual effective mass approximation. However, for 
nanocrystals, the effective mass itself becomes size dependent and the exponent X has 
been reported to vary from 1.2 - 1.8 and becomes 2 at large size^^^'^^^. For the purpose of 
our simulation we consider two cases 

a) X = 2, C = 486 eV/A^ (Ego = 1-71 eV for d = 30 A) (22,24) _ 
and 

b) X = 1.4, C = 126 eV/A^-^ (Ego - 2.25 eV for d = 30 A) (2^) . 

The interaction volume of the column structures with incident light is V(d) — 
Ld^ . The resulting distribution R(d) = P(d) x V(d) is properly normahzed to have equal 
number of dipole oscillators in each case. We now make a change of variable from d to Eg 
( using Eq = Eg + C/d^) to get the corresponding distribution P(Eg) of the band gaps. 

For a Gaussian distribution of sizes (eq.9) - the distribution of band gaps is 

P'^iEa) = PoiEa - ^j(-^-3)Aea;p[-A((|^^^)V^ - 1)^] (12) 

where Pq is a normalization constant, A — 0.5((io/cr)^ and E^ is the maximum value of Eg. 
Following a similar procedure for the Lognormal distribution (eq.lO) we have 

P^(Ea) = MEa - E,)<--.«.r'e.p[-B(^|^g^ - 1)^] (13) 

where Rq is the normalization constant , B = 0.5 (mo/iTL)^ . 

Fig. 4 shows the distribution of band gaps P(Eg) when the energy upshift is 
given for X = 1.4 and X = 2. For the Lognormal case, we see that P(Eg) peaks close to 
the band gap of c-si. 

Equations 8, 12 and 13 are used to calculate a{E) for each incident energy E = 

hu!. 

a{E) for a Gaussian distribution of sizes is given by 
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(E) = aolJ\E - Ea)-^{Ea - E,)^--^-'^^ exp[-A{{^-^)y^) - ir]dEa (14) 



where ao is a normalization constant . 

For a Lognormal distribution of sizes a{E) is 



(15) 

where /3o is a normahzation constant. 

Equations (14) and (15) are used to simulate various possible absorption 

processes . 



IV. RESULTS OF THE SIMULATION 

In this section, we present the results for the absorption co-efficient as a 
function of energy for different cases on the basis of the above model. The values of various 
exponents for different physical situations are summarized in Table I. 
(1) Direct band gap 

Figure 5 shows a{hu!) vs hcu for a direct band gap material where k is conserved 
in the optical transition and no size dependence of the oscillator strengths ( 7 = ) is 
considered (cases 1 and 3 and cases 5 and 7 in Table 1 ). To facihtate comparison between 
various cases , we have normalized the absorption to its peak value in each case. We see 
from the figure that the absorption tends to fall off at high energy ( in sharp contrast with 
the experimental data ). These results of the simulation can be easily understood as 
follows. 

For a ID system, the JDOS has a singularity at E = Eq and falls off for E > 
Eq. The dominant contribution (to each ck{E)) is hence only from a particular set of 
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nanostructures with band gap Eq where E —Eq and all other Eqs have almost no 
contribution. So the simulated a{E) vs E plot mimics P(-E'g)- We hence conclude from the 
simulation and the experimental data that PS is not a direct band gap material and we 
will not consider this case further. 
(2) Indirect band gap 

We now consider the case when momentum(k) is not conserved for the optical 
transition. Here we discuss two cases - size independent oscillator strengths (7 = 0) and 
size dependent oscillator strengths ( assuming 7 = 6). 

We first consider the situation where the oscillator strength has no size 
dependence. Figure 6 (curves (a) and (c)) is a plot of the calculated a vs hu; for a Gaussian 
distribution of sizes (case 2 and 4 of table I ). Figure 7 (curves (a) and (c)) is a plot of the 
calculated a vs huj for a Lognormal distribution of sizes (case 6 and 8 of Table I). Also 
shown in the figures is a plot of the experimental curve for a p type PS. (The dip in a near 
2.2 eV in the experimental curve is an artifact related to a filter change). To facihtate 
comparison , the calculated a is scaled to match the experimental curve near 3 eV. We see 
from the plots 6(a) and 6(c), that a initially increases with energy and then saturates - at 
variance with the experimental results. We understand this result qualitatively as follows. 
In section 1 , we have already mentioned that the JDOS is a constant independent of 
energy for hcu < Eq and zero for huj > Eq- Hence for hcu > Eq ,a given nanostructure 
contributes equally at all energies to a. However P(Ec;) falls off as E^ increases. 
Nanostructures having a band gap {Eq) above a certain cut off in energy (in the high 
energy tail of the distribution) do not contribute significantly to a{l%uj) (as P(Eg) is very 
small). This gives rise to the saturation of a{E) in figures 6(a) and 6(c). This saturation is 
more obvious in figures 7(a) and 7(c), because P(Eg) falls off very rapidly for a Lognormal 
distribution of sizes. These results also do not agree with the experimental data. 

The situation is not remedied even if we include higher excited states (n-^ = 2, 
Tiy = 1\ etc.) in the energy range considered in our calculation. Based on the above 
discussion, we see that for a to increase with fiuj, it is necessary to increase the contribution 
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to a from the nanostructures in the high energy region of a{E) .This is possible if the 
oscillator strength is enhanced for the nanostructures with large Eq ( smaller sizes ) . 

We now consider the case when ¥[Eg) varies as f /(P . Figure 6 (curves (b) and 
(d) for case 2 and 4 in Table I )and figure 7 (curves (b) and (d) for case 6 and 8 Table I ) 
are plots of a{huj) vs Tiui for 7 = 6. For the Gaussian distribution that we have chosen, we 
see that the absorption still saturates although the onset of saturation moves to slightly 
higher energy. This is a measure of the contribution of the (high energy) tail of P{Eg) for 
the Gaussian case to the total absorption. For the Lognormal distibution, we see from 
figure 7 (curves (b) and (d)) that a{liuj) vs Uuj now qualitatively looks like the 
experimental curve. 

From a comparison of figures 6 and 7, we find that in order to get the 
absorption co-efiicient to increase with incident energy , it is important that the 
contributions to a from the low energy regions of P(£'g) do not dominate the optical 
absorption at large hu. This is also satisfied if the distribution of band gaps ^{Eg) is 
sufficiently broad. We illustrate this in Fig. 8 where we plot the calculated value for a for a 
Gaussian distribution of band gaps with the mean energy at 2.25 eV and a variance of 0.75 
eV. Hence we find that the nature of the band gap distribution P{Eg) plays a crucial role 
in determining the shape of the absorption spectrum . 

We conclude that to properly account for the optical absorption in PS we need 

to assume 

(1) PS is a pseudo ID indirect band gap semiconductor. 

(2) The oscillator strength is size dependent and is given by a power law f /(P . 

(3) The distribution P{Eg) also plays an important role . Lognormal 
distribution of the nanostructures diameters in porous silicon seems to naturally account 
for the absorption spectrum. 

All these considerations are in general agreement with results of luminescence 
experiments on PS .The low temperature luminescence of PS excited near the peak of the 
broad luminescence spectrum exhibits steps which have been associated with the 
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characteristic TO mode of c-silicon . More recently*^^°) well resolved peaks 
corresponding to the TO mode of Si have been reported for luminescence of PS suggesting 
that PS is an indirect band gap semiconductor. The observed dominant zero phonon 
component of luminescence constitutes direct evidence for the enhancement of the 
oscillator strength due to strong quantum confinement at smaller sizes . 

We now comment on the validity of the procedure for obtaining the band gap of 
PS from \J olTuij vs hio. Figure 9 is a simulated plot of \J ahio vs hio ( case 8 in Table I ) for 
three different values of 7 and the experimental data on p type sample . First, we note 
that the intercept of the linear portion of each of the simulated curves is neither coincident 
with the peak position ( near 1.3 eV ) of the corresponding band gap distribution ( see 
figure 4 ) nor close to the mean band gap energy £'0 — 2.2b ev. So we see that it is very 
difficult to relate this intercept with such physical quantities .Therefore , it is not obvious 
how to interpret the intercept with some kind 'effective band gap '. We also note that the 
intercept depends on the choice of 7. We hence believe that the intercept does not have 
the simple interpretation of a hand gap as defined for a homogeneous system. We would 
like to point out that for such a heterogeneous system the gap can be interpreted only in 
an operational sense like E03 or E04 - depending on the value of the absorption coefficient 
(after correcting for the sample porosity). 

Finally we try to understand the temperature dependence of the band gap. To 
explain the rigid shift of the absorption curve with temperature (to lower a), we argue 
that Eq of each of the nanostructures follows a relation Eq (T) = Eqq — ■jT, where Eqq is 
the band gap at zero degree Kelvin . We assume that the change of the band gap with 
temperature is a consequence, like in bulk semiconductors, of the electron -phonon 
interaction and take 7 = 0.5 meV/K. In our model, changing temperature is equivalent to 
changing the value of C in the expression that determines the energy upshift AE. Figure 
10 is a plot of ln(Q;) vs hoj for three different values of C. We see that the absorption curves 
are rigidly shifted vertically with respect to each other as also seen in figure 2. We 
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argue that this provides a natural explanation for the observed temperature dependence of 
a . 

V. CONCLUSIONS 

In this paper we have investigated and clarified the absorption process in free 
standing thin films of PS . 

We have elaborately discussed that one can not use eq 1. to obtain the band 
gap of porous sihcon . We point out that these objections are also true for other low- 
dimensional nanostructures like quantum wires , quantum dots etc . One can easily 
generalize our results ( eq 8 , 14 , 15 ) and use them to analyse the absorption process for 
other nanostructured materials too . For that purpose one just have to put appropriate 
density of states function and the size distribution in those equations . 

We have calculated the absorption coefiicient as a function of energy assuming 
that PS can be modelled as a pseudo ID system having a distribution of band gaps . 
Prom our calculations, we show that it is neccessary to look at PS as an indirect band gap 
material. To satisfactorily account for the absorption, we need to invoke that the 
oscillator strengths has a size dependence. We also show that a Lognormal distribution 
of diameters of the nanostructures that make up PS appears to account for the measured 
absorption spectrum of p PS . 

We can explain the temperature dependence of the absorption in PS on the 
basis of our model too. 

We have also shown that porosity can be inferred non — destructively from a 
transmission measurement in the region of low absorption . 
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Gaussian distribution in d 


Lognormal distribution in d 


Case 


Value of X 


Nature of Transition 


Case 


Value of X 


Nature of Transition 


1. 


2.0 


Direct 


5. 


2.0 


Direct 


2. 


2.0 


Indirect 


6. 


2.0 


Indirect 


3. 


1.4 


Direct 


7. 


1.4 


Direct 


4. 


1.4 


Indirect 


8. 


1.4 


Indirect 



TABLE I. These physical situations we have considered in the model for the optical 
absorption in porous silicon .For direct transitions W = 0.5 and for indirect transitions W 
— . The value of 7 is taken as 6.0 for size dependent oscillator strengths and zero for the 
size independent case. 
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FIGURES 

FIG. 1. Porosity calculated ( using eq. 7 ) from the transmission curves*^^'-*^^) are plotted 
against the corresponding gravimetrically determined values of porosity . The points fall nearly on 
a straight line of slope one. (a) data from ref. 9 and (b) data from ref. 16 . 

FIG. 2. Temperature dependence of the absorption spectrum of p type porous silicon . (a) 
300 K and (b) 100 K . 

FIG. 3. Size distribution of the column widths . (a) For a Gaussian and (b) for a Lognormal 
distribution of widths, having do=30 Aand a= 4 A. 

FIG. 4. Band gap distribution for corresponding size distribution given in figure 4. (a) and 
(b) are for Gaussian distribution for X = 2 and 1.4 respectively, (c) and (d) are for the Lognormal 
distribution for X = 2 and 1.4 respectively. 

FIG. 5. Simulated a{E) vs E plot for direct transitions (7 = 0). Plots (a) and (b) are for 

the cases (1) and (3) in Table I . Plots (c) and (d) are for the cases (5) and (7) in Table I. 

FIG. 6. Simulated a{E) vs E plot for indirect transitions .(a) and (b) correspond to case (2) 
in table I for 7 = and 6 respectively .(c) and (d) correspond to case (4) in table I for 7=0 and 6 
respectively .(e) Measured a{E) vs E plot for p type PS. 

FIG. 7. Simulated a{E) vs E plot for indirect transitions .(a) and (b) correspond to case (6) 
in table I for 7=0 and 6 respectively .(c) and (d) correspond to case (8) in table I for 7=0 and 6 
respectively . (e) Measured a{E) vs E plot for p type PS. 

FIG. 8. (a)Simulated a{E) vs E plot for broad gaussian distribution of band gaps ( Eq = 2.25 
eV and cte = 0.75 eV ) with indirect transitions , X=1.4 and 7 = 6. (b) Measured a{E) vs E 
plot for p type PS. 

FIG. 9. Simulated plots of (a?iw)^/^ vs hcv for (a) 7 = 5 , (b) 7 = 5.5 , (c) 7 = 6.0 and (d) 
corresponding plot for the experimental data on p type PS . 
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FIG. 10. Simulated log(a(£')) vs E plots for various temperatures ( for case 8 with 7 = 6). 
(a) C = 126 = Co , (b) C = 1.25Co , (c) C = 1.5Co . Curves are rigidly shifted vertically from one 
another. 
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